Resistive (/-modes in a reversed field pinch plasma 
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First direct experimental evidence of high frequency, high toroidal mode number (n > 20), mag- 
netic fluctuations due to unstable resistive interchange modes (g-modes) resonant in the edge 
region of a reversed fleld pinch (RFP) plasma is presented. Experimental characterization of 
time and space periodicities of the modes is provided by means of highly resolved in-vessel edge 
and insertable magnetic diagnostics. It is found that the spectral mode properties are in good 
agreement with the predictions of the theoretical linear resistive magnetohydrodynamic stabil- 
ity analysis. A simple model is proposed for the observed saturation levels of the modes. Pacs: 
52.25.Xz,52.30.Cv,52.35.Py,52.55.Hc,52.65.Kj,52.70.Ds 



MHD turbulence of interchange nature is ubiquitous in 
both astrophysical and laboratory plasmas. In fusion ori- 
ented magnetically confined plasmas interchange modes 
(also known as g-modes) can cause significant anomalous 
and particle transport across the magnetic field. In toka- 
mak plasmas g-modes are counteracted by shaping and 
curvature effects of the mean field. The remaining tur- 
bulence is therefore of ballooning nature being espe- 
cially concentrated in the unfavorable toroidal curvature 
region. For example. Edge Localized Modes (ELMs) are 
interpreted in this way [^. In the reversed field pinch 
(RFP) 13 the average edge curvature is unfavorable, be- 
ing the field lines mainly in the poloidal direction, there- 
fore interchange modes should be expected. According to 
Suydam's criterion [3| the effect of unfavorable curvature 
can be mitigated by enhancing the magnetic shear, and 
this is particularly relevant for RFPs. However, finite 
plasma resistivity can give rise to a resistive interchange 
branch of instabilities, known as resistive g-modes, that 
can be destabilized well below the Suydam's limit. In the 
existing literature (/-modes have been often envisaged to 
be linearly unstable and invoked as responsible for mag- 
netic field ergodization and enhanced diffusion and trans- 
port in RFP plasmas 0, IH 0, [!]■ Some experimental 
results on this subject are described in ref. In this 
Letter, we present the first clear experimental evidence 
of the existence within the magnetic spectrum of resistive 
g-modes characterised by high mode numbers along with 
a first direct comparison with modeling prediction. 

The experimental activity here described has been per- 
formed on RFX- mod the largest operating toroidal 
RFP device (major radius i? = 2 m, minor radius 
a = 0.459 m). Hydrogen is the working gas and the 
plasma current Ip spans between 0.3 and 1.8 MA. The 
on-axis electron temperature at the highest plasma cur- 
rent exceeds 1.5 keV, while the volume averaged electron 



density Ue is in the range 1 — 10 x 10 
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duration is up to 0.5 s. The analysis is performed on 
a large number of discharges (~1000) with reversal pa- 
rameter, defined as F = B^{a)/ < > {B^{a) and 
< B^ > toroidal field at the wall and toroidal field aver- 
aged over the plasma cross-section, respectively), in the 
range [—0.5,0], which corresponds to edge safety factor 



q{a) G [—0.07,0]. The plasma is ohmically heated and 
the loop voltage is in the range 20 - 40 V, which gives to 
the RFP the attractive feature of being able to use ohmic 
heating only also in fusion-relevant conditions. In the 
RFX-mod device the magnetic boundary is determined 
by a thin Cu shell, with vertical field penetration time 
of 50 ms, located at r/a = 1.12, and by a system of 192 
active saddle coils, supervised by a digital feedback sys- 
tem, fully covering the machine with the aim of control- 
ling the radial fields due to field errors and MHD modes 
[ll| . In the experiments described herein the mesh of 
saddle coils has been used in the so-dubbed clean-mode- 
control (CMC) configuration, in which the flux through 
each sensor loop is controlled by the corresponding active 
coil after performing the two-dimensional EFT of Br and 
Bcf, measurements and computing a real time correction 
of the aliasing of the sideband harmonics generated by 
the discrete saddle coils [l2j . 

Magnetic fluctuations have been measured by two dif- 
ferent highly space and time resolved systems of in- 
vessel magnetic probes: the first is an insertable edge 
probe, used to investigate high order toroidal harmonics 
(|7i| < 85), the second consists of arrays of probes, cover- 
ing the full toroidal and poloidal angles of the torus. The 
insertable probe, named U-probe [IJ], which was used in a 
set of experiments performed at low toroidal plasma cur- 
rent {Ip < 500 kA), consists of two Boron Nitride cases, 
5 cm toroidally spaced, radially inserted from the vac- 
uum chamber up to r/a « 0.9 at = 0° (defined on the 
equatorial plane on the low field side of the machine with 
positive 9 pointing in the upward direction), without sig- 
nificant perturbation of the plasma. Each case contains, 
along with a number of electrostatic pins, a radial array 
of 7 three- axial magnetic pick-up coils measuring the time 
derivative of the three components of the magnetic field 
{Br, Bg, B^). The sampling frequency is 10 MHz, with 
an estimated bandwidth up to 3 MHz. The second sys- 
tem of probes is a subset of the ISIS (Integrated System of 
Internal Sensors) diagnostics [13], a large set of high fre- 
quency electrostatic and magnetic probes located inside 
the vacuum vessel. The magnetic probes consist of coils 
measuring B^j,, placed at r/a = 1.03 behind the graphite 
tiles which cover the first wall of the machine. The probes 
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FIG. 1: Typical power spectra of the time derivative of the 
three components of the magnetic field, measured by one of 
the three-axial coils housed in the U-probe. 



used in this study are distributed in two different arrays: 
the first constituted of 48 coils evenly distributed in the 
toroidal direction located a,t 9 — 250°; the second one 
is a poloidal array, made of 8 equally spaced coils. The 
sampling frequency for the magnetic measurements from 
the ISIS system is 2 MHz, while the estimated bandwidth 
is up to 400 kHz. The two arrays allow a resolution of 
the toroidal and poloidal mode numbers n and m up to 
24 and 4, respectively. 

In Fig. [T]the typical power spectra of the signals from 
the three components of the magnetic field fluctuations 
{Br, Bg, B^), taken with the U-probe in a deeply re- 
versed {q{a) « —0.04) discharge, are shown. The spectra 
show that the highest fluctuation levels are present in the 
(Br,B(j,) signals, which correspond to the perpendicular 
components at the edge of the RFP configuration, where 
the dominant magnetic field is Bg . The spectra show that 
large magnetic fluctuations are present for frequencies 
up to about 200 kHz, with a large part of the magnetic 
fluctuating energy concentrated in the frequency region 
spanning from 30 to 150 kHz. The space-time proper- 
ties of the magnetic fluctuations have been obtained by 
Fourier decomposing the signals coming from the com- 
plete toroidal array of coils measuring i?^. In Fig. [2] the 
frequency spectrum of each toroidal Fourier component 
is plotted in a color-coded plot for two different exper- 
imental operating conditions, corresponding to plasma 
equilibria characterised by the q profiles plotted in Fig. 
Ht. The q profiles are deduced by using a suitable mag- 
netic profile reconstruction (see for example [T^). Each 
of the spectra has been obtained as an average over about 
200 discharges. Fig. [2^ has been obtained by analysing 
data taken in a set of discharges with deep reversal of the 
toroidal magnetic field (g(a) ~ —0.05); the S{n,f) spec- 
trum exhibits the excitation of coherent magnetic fluctu- 
ations, in the form of a broad peak, at frequencies above 
30 kHz and high (n up to 50) positive toroidal mode num- 
bers (in the RFX-mod convention, negative and positive 
n values refer to modes resonant inside and outside the 




FIG. 2: Typical color coded contour plot of the S{n, f) spec- 
tra of magnetic fluctuations for two different experimental 
conditions: a) S{n, f) spectrum obtained at q{a) « —0.05; b) 
S{n, /) spectrum obtained at shallower reversal of the mag- 
netic field {q{a) ~ —0.015); c) the two reference q{r) profiles 
(the vertical blue line marks the plasma edge). In the sets of 
discharges considered 300kA < Ip < GOOkA. 



reversal surface, respectively). It is worth noting that due 
to the limited (l^il < 24) spatial resolution of the toroidal 
array of coils, the computed S(n, f) spectra, originally 
affected by aliasing, had to be corrected. This opera- 
tion was made possible by a comparison to the S(n, /) 
spectra deduced by applying the two-point technique [l8| 
to the signals from the closely spaced magnetic coils of 
the U-probe, which show that the observed peak actually 
correspond to high positive n values. Fig. [Jb shows the 
result of the same analysis performed on discharges char- 
acterised by a shallow reversal of the toroidal magnetic 
field, q{a) « —0.015. In this case no peak is present and 
the magnetic fluctuations are concentrated at frequencies 
below 50 kHz, positioned on an almost continuous linear 
dispersion relation, in the form lo — -^Vph, with uj = 27r/, 
where / is the frequency, n is the toroidal mode num- 
ber, and Vph is the phase velocity of the fluctuations in 
the toroidal direction. The measured phase velocity is of 
the order of 20-30 km/s in the counter toroidal plasma 
current direction, which is in good agreement with the 
plasma flow velocity at the edge of the RFX-mod device 
deduced by means of a gas puffing imaging diagnostics 
[l^ and is also consistent with the E x B flow measured 
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FIG. 3: a) Theoretical growth rate for different toroidal mode 
numbers, corresponding to two q{a). Different symbols mark 
the two different plasma pressure profiles used for the simu- 
lation, b) Experimental S{n) spectra in the same q{a) condi- 
tions of a) (experimental spectra are obtained as average over 
10 comparable discharges, in terms of Ip, rie and q{a)). 



by Langmuir probes . It is important to note that the 
peak in the {n, f) plane of Fig. [^h- results to be almost 
aligned with the linear dispersion relation characterising 
the continuum low frequency part of the spectrum. This 
is an indication that the frequency associated to the mode 
with n > 20, is mainly due to a Doppler effect, as the 
mode rotates with the plasma. The analysis for different 
values of q{a) reveals that the value of n corresponding 
to the maximum of the fluctuation amplitude is observed 
to follow the relation n ■ \q{a)\ «1. The analysis of the 
poloidal periodicities of the magnetic signals from the 
poloidal array of coils reveals a clear m — 1 nature of the 
magnetic fluctuations in the explored frequency range. 
The measured poloidal and toroidal mode numbers thus 
indicate that the observed high frequency modes corre- 
spond to magnetic perturbations resonating at the edge 
of the plasma column, in the region between the reversal 
surface and the inner surface of the graphite tiles con- 
stituting the first wall of the RFX-mod machine. In the 
RFP nomenclature these are called externally resonant 
modes (ist - 

The measured spectra have been compared to those 
predicted by a linear stability analysis performed by 
means of the cylindrical code ETAW, already extensively 
used for RFP calculations fl5| and more recently suc- 
cessfully benchmarked against the MARS code . The 
code solves the linear cylindrical resistive incompressible 
and inviscid single fluid MHD equations, using a spectral 
formulation and a matrix shooting eigenvalue scheme. 
Up to two resistive walls are considered for the bound- 
ary conditions, with a thin shell approximation. The 
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FIG. 4: Total normalised amplitude of the magnetic pertur- 
bation as a function of the Lundquist number S for a given 
plasma equilibrium {q{a) = —0.03). Dashed line, referring to 
the right hand side y-axis, represents the dependence on 5* of 
the growth rate predicted by the ETAW code for the same 
plasma equilibrium. 



plasma model is solved inside the flrst wall and the so- 
lution is then matched to the external solution of the 
vacuum cylindrical Laplace equation, analytically known 
in terms of Modifled Bessel Functions. For the problem 
under consideration we have assumed that only one wall 
is present, which corresponds to a perfect ideal shell at a 
shell proximity (b/a) of about 1.05, i.e. 5 % of the plasma 
minor radius. This choice appears to be the most appro- 
priate, since, as it has been shown by the measurements, 
the modes rotate at relatively high speed together with 
the plasma, and, on the other hand, a metal liner with 
1-2 ms magnetic field penetration time, is present in the 
RFX-mod device. Therefore the liner acts as a perfect 
conductor for these fast rotating modes. 

With these assumptions, we obtained the spectrum 
of unstable modes given in Fig. [3^ for two different 
q profiles, one with a deeper reversal, corresponding to 
q{a) = —0.05, in which the spectrum is peaked at a lower 
n value (n « 25), and one at a shallower reversal, with 
q{a) = —0.03 , with a spectrum peak at n w 40. The 
/3 used in the simulation is 2%, comparable to the ex- 
perimental one for the considered case. The theoretical 
predictions are in good agreement with the experimental 
toroidal mode number spectra S{n), shown in Fig. [3]d, 
obtained by integrating the S{n, /) spectra over frequen- 
cies above 30 kHz, which, as discussed above, are those 
pertaining to the investigated m = 1 modes (in order to 
obtain the value of the fluctuating magnetic fleld, mag- 
netic signals have been numerically integrated). In par- 
ticular, similarly to the experimental results, the toroidal 
mode number spectrum for a given equilibrium of Fig. [3^ 
shows that the growth rate 7 (normalised to the Alfven 
time) vanishes for modes resonant close to the stabilising 
wall (at low n's), and also for modes resonant close to the 
reversal of the toroidal magnetic fleld (at high n's), where 
the effect of the magnetic shear is stronger. As also shown 
in Fig. [3^, the theoretical growth rates at a given f3 are 
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FIG. 5: Total (normalised) magnetic perturbation 5B^ in- 
duced by the resistive g-modes as a function of q{a) (full cir- 
cles). Theoretical growth rate, normalized to the squared n 
values as a function of q{a) (open circles). 



strongly influenced by the pressure profile. In particular, 
larger edge pressure gradients are associated to larger 7 
values. Moreover, modes are found to be linearly sta- 
ble at zero f3. The model also indicates that the growth 
rates scale like ~ S~^'^ {S being the Lundquist number), 
as shown in Fig. |31 for a given plasma equilibrium cor- 
responding to q{a) — —0.03, which is very close to the 
theoretically expected S"^^/"^. In the same figure, the the- 
oretical growth rates are compared to the experimental 
total amplitude of the magnetic fluctuations, obtained by 
summing the contribution from each single toroidal har- 
monic. Despite the large dispersion of the experimental 
points, a scaling close to S*^^/"^ is fully compatible with 
the data. 

For all these reasons, i.e. effect of the resonance, ef- 
fect of the beta and shear and Lundquist scaling, we can 
classify these modes as resistive g. These are interchange- 
like instabilities which can develop in a RFP at relatively 
low P values [l^l , well below the Suydam's limit for ideal 



interchanges. 

As a last point we want to address the problem of 
the nonlinear saturation of these modes. Clearly for 
this purpose a nonlinear model would be required (see 
for example [2l|). and/or a model taking into account 
all the effects which could affect the stability proper- 
ties of g-modes extensively described in the literature, 
with the inclusion of finite Larmor effects, parallel ion 
viscosity, diamagnetics effects [H, HI] (and references 
therein). However, we consider a simple model in which 
the growth rate is balanced by a dissipation mechanism 
proportional to the squared mode wave-number, which 
is the consequence of having a Laplacian-like dissipation 
term in the equations. In this case, it is expected that 
the measured fluctuation amplitude correlates with some 
quantity which takes into account both the growth rate 
and the mode number. As shown in Fig. [51 despite the 
semplicity of the model, a very nice correlation is found 
between the maximum growth rate divided by the square 
of the mode number (corresponding to the maximum) at 
different q{a) values and the experimental amplitude of 
the fluctuations. In particular, in Fig. Othe amplitude 
of the total perturbation produced by these modes is ob- 
served to reach values up to almost 1% of the equilibrium 
magnetic field at the edge in deeply reversed discharges. 
It is worth to note that the magnetic perturbation pro- 
duced is only a factor four smaller than that due to the 
core resonant dynamo modes, which are essential to sus- 
tain the RFP configuration [l0| . Therefore we think that 
the described phenomenon could play an important role 
in determining the edge transport properties in RFP plas- 
mas, for which the exact mechanism is still under debate. 
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